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We consider the one-dimensional half-filled Hubbard model. We show that 
the excitation spectrum is given by the scattering states of four elementary 
excitations, which form the fundamental representation of SU{2) x SU{2). 
We determine the exact two-particle Scattering matrix, which a solution of 
the Yang-Baxter equation and reflects the <S'0(4) symmetry of the model. 
The results for repulsive and attractive Hubbard model are related by an 
interchange of spin and charge degrees of freedom. 
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1. Introduction 

Strongly correlated electronic systems are currently under intense study in relation 
with the phenomenon of high-T^ superconductivity. The two-dimensional Hubbard model is 
believed to be the most promising candidate for an electronic theory of superconductivity. 
It is believed to share important features with its one-dimensional analog'^l. An important 
problem in understanding the dynamics of the Hubbard model is the separation of spin and 
charge in one'^'^l and two dimensions^ . A particularly interesting feature of the Hubbard 
model is that in one dimension it can be solved exactly by means of the Bethe Ansatzl^l . Thus 
it is possible to obtain a complete and unambiguous picture of the dynamics in one dimension. 
The main purpose of this paper is to use the Bethe Ansatz solution to demonstrate that the 
one-dimensional Hubbard model is a factorizable scattering theory of four quasiparticles. 
These elementary excitations transform in the fundamental representation of SU (2) x SU (2). 

The Hubbard hamiltonian is given by the following expression 

^ = - E E (sU+i,. + + 4t^EKT - - \) ■ (1-1) H 

Here c] ^ are canonical fermionic creation operators on the lattice, j labels the sites of a 
chain of even length L, a labels the spin degrees of freedom, U is the coupling constant, and 

= '^j.o'S.o- number operator for spin a on site j. The hamiltonian (1.1) is invariant 

under a 'so'(4)=SU(2)xSU(2)/Z2 algebra'^-^^l The two SU(2) algebras are the spin SU(2) 
generated by 

L L L ^ 

= E sItSu ' = E 4iS-,T ' = E 2 - ".-.t) ' (1-2) spin 

j=i j=i j=i 

and the ?7-SU(2) algebra 

^ = E(-l)S-,TSa' = E(-l)Ma4T ' ^^ = ^EKT+^.i-l)- (1-3) eta 
j=i j=i j=i 

All 6 generators commute with the hamiltonian (1.1). Clearly the operator + r]^ has only 
integer eigenvalues (as L is even) and all half-odd integer representations of SU (2) x SU (2) 
are projected out. Thus the symmetry algebra is 5*0(4) as asserted above. 

The nested Bethe Ansatz for the Hubbard model'^' provides eigenstates of the hamil- 
tonian that are parametrized by sets of so-called spectral parameters kj and A^. These 
parameters are subject to the Lieb-Wu equations 

.. r -TT sinik.:) — K„ + iU 
sinikA - A„-iU 

(1-4) liebwu 

sinjk,) -K, + iU _ A, - A, + 2iU 

11 sin{k,)-K,-iU l}^K-K,-2iU ' ^ '-' ' 
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These equations are the basis for the determination of ground state and excitation spectrum 
in the thermodynamic hmit. 

There has been substantial previous work on the excitation spectrum of the repulsive, half- 
filled Hubbard model. A. A. Ovchinnikov calculated the spin-triplet excitation'^^1, which was 
then re-examined by T.C Choy and W. Young in [12]. In his series of excellent papers'^'^'^^1 
F. Woynarovich gave a very detailed analysis of both spin and charge excitations. The at- 
tractive case was studied by B. Sutherland, who investigated the Bethe Ansatz structure 
of the ground state, and by F. Woynarovich, who determined ground state and excitation 
spectrum from the known results for the repulsive case by using discrete symmetries of the 
hamiltoniant^^l . This method, while providing a nice picture of the relations between the 
Bethe Ansatze for repulsion and attraction has the drawback of making assumptions about 
transformation properties of spectral parameters under the map connecting the [/ > and 
[/ < cases. In section 2 we present a direct Bethe Ansatz analysis of the ground state and 
excitation spectrum in the two-particle sector of the attractive Hubbard model. Our direct 
computation as expected confirms Woynarovich's results for the dispersion relations'^^' . The 
main purpose of this section is to determine the so-called "shift-functions" I^^'^^l , which are 
needed to evaluate the two-particle S-matrix. In section 3 we give a quasiparticle inter- 
pretation of the excitation spectrum. All excited states are scattering states of only four 
quasiparticles (or elementary excitations), which form the fundamental representation of 
SU{2) X SU{2). Two of these elementary excitations carry spin but no charge, and two 
carry charge but no spin. This reflects the separation of spin and charge degrees of freedom 
in the one-dimensional Hubbard model. In section 4 we evaluate the exact S-matrix of the 
attractive Hubbard model. It is a solution of the Yang-Baxter equation, which implies the 
two-particle reducibility of the A^-body S-matrix. In section 5 we give the results for the 
repulsive case, and in section 6 we show that attractive and repulsive Hubbard model are 
related by an interchange of spin-and charge degrees of freedom. This relation holds both on 
the level of the quasiparticle spectrum and on the level of the S-matrix. Finally section 7 is 
devoted to a summary and discussion of our results. The appendix deals with the problem of 
showing that all excitations in the A^-particle sector can be interpreted as scattering states 
of the four elementary excitations. 



2. Attractive Hubbard Model 

In this section we construct the ground state and all excitations in the two-particle 
sector for the attractive case U < 0. Starting point are the Lieb-Wu equations (1.4). The 
solutions of (1.4) can be split into three different types of subsets (strings), which are^-'^'^l 

(1) a single real momentum 

(2) m A^'s combine into a string-type configuration ('A-strings'); this includes the case 
m — 1, which is just a single real A^, 
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(3) 2m A;j's and m Aq,'s combine into a different string-type configuration ('A;-A-strings') 
For a A-string of length m the rapidities involved are 



K'' = A^ + i(m + 1 - 2j) I [/| A^ real j = 1, 2, . 
The A;'s and A's involved in a A;-A-string are 

£ = sin-^(A;™ + zm|[/|) 

sin-HAr + ^(m-2)|t/|) , = tt - /c^ 
sin-^(Ar + i(m-4)|t/|) , = tt - fc^ 



and 



k 



2m-2 



.2m 



sin-^(A;^'"-i(m-2)|C/|) , kl""-^ 
sm-'iA'^"" - im\U\) . 



TT — k 



2m-2 



(2.1) lambdastr 



(2.2) klstrl 



A'^^'^ ^A'^^ + i{m + l-2j)\U\ , A'^"*real j = l,2,...,m. 



(2.3) klstr2 



Equations (2.1) - (2.3) are valid up to exponential corrections of order C(exp(— 5L)), where 
5 is real and positive (as long as the real parts of the spectral parameters are much smaller 
than L) . Note that the /c-A-strings for the attractive Hubbard model are different from the 
repulsive case treated by M. Takahashi in [17]. 

We now consider a solution of (1.4) with M„ strings of type (2.1) of length n, M'^^ k-A- 
strings (2.2), (2.3) of length n, and a total number of electrons A^g. Inserting the prescription 
(2.1) - (2.3) into the Lieb-Wu equations (1.4) and taking the logarithm we arrive at 



11=1 a=l ^ 11/ n=l a=l 



sin kj — A'^ 



N^-2M' . , 

/ - sm kj 
U V n\U\ 



CO M„ 

^^Ja + E E 
m=l 13=1 



n\U\ 



An _ Am 



L [sin-^ {A!^ + in\U\) + sm"^ (A;," - in\U\)\ 

N^-2M' , . , 

sm k. 



ML 



2-^?+ E »(^^9?^) + EEe™. 

I I / m=l /3=1 



A /n A /m 



(2.4) pbc 

(2.5) pbclb 

(2.6) pbc2 



where L — 2 x odd is the even length of the lattice, J , J", and J'^ are integer or 



half-odd integer numbers, M' — S^i^-^n; ^(^) — 2 arctan(a:), and 



^ |n— ml / ^|n— m|+2/ V 



|n-m|+2^ 



^ n+m—2 ■ 



X 



n+m 



2^(f) + 2^(f) + ... + 2^(^) + ^(|-) 



) if n 7^ m 
if n = m 



(2.7) theta 
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The integer (half-odd integer) numbers in (2.4)-(2.6) have the ranges 

^ oo 
m=l 

^ oo 

\ J'a \<7;(L - + 2M' - J2 *„™< - 1) . 



< L < 



2 - ^ - 2 

m=l 



where = 2min{m, n} — S^^. Each set of "integers" {Ij}, {Ja}^ {^a} in one-to-one 
correspondence with a set of spectral parameters, which in turn unambiguously specifies 
one eigcnstatc of the hamiltonian (1.1). Energy and momentum of a solution of the system 
(2.4)-(2.6) are found to be 



N^-2M' oo 

cos(/t;) -4^^Re|Vl- (Aa" + «|^|?^)^} -2f/iV + L/L , (2.9) E 

1=1 n=l a=l 

N^-2M' oo 

j=l n=l a=l 

Equations (2.4) - (2.10) can now be used to determine the ground state energy and to con- 
struct excitations. The ground state can be found by taking the zero temperature limit of the 
Thermodynamic Bethe Ansatz based on (2.4)-(2.9), which can be constructed analogously 
to the repulsive case^^^^ (see the Appendix of this paper). An important element in our 
construction of excitations over the ground state is a theorem proved in [18], which states 
that 

(i) All eigenstates of the hamiltonian (1.1) with finite spectral parameters, < and 
N^< L are lowest weight states of the 5*0(4) algebra (1.2), (1.3). These states are called 

regular Bcthc states. 

(ii) The set of states obtained by acting with the SO{A) raising-operators on the regular 
Bethe states is complete and forms a basis of the electronic Hilbert space. 

This theorem implies that excitations constructed by means of the Bethe Ansatz arc all 
lowest weight states of the SO (4) algebra. Furthermore all S'0(4)-descendants of these 
excitations must be taken into account in order to obtain a complete set of excited 
states. 

Ground State 

The ground state for the half- filled band is characterized by choosing M[ — ^, = 
Vn, = Vn > 2, and fiUing all ^ vacancies given by (2.8) for the integers J'^ sym- 
metrically between and The Lieb-Wu equations (2.4)-(2.6) for the ground state 
are 

2LRe{arcsin(A, + ^|t/|)} = 27rj;^ + ^^(^^|^) , a=l,...,^. (2.11) GSEQ 

p=i I I 
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Subtracting equations (2.11) for consecutive a's and using J^^i ~ J'a — ^ we obtain an 
equation for the finite-interval density 

Po(Aa) = T(x ^ X ^ ■ (2-12) rho 

In the thermodynamic hmit this equation turns into an integral equation for the ground 
state density Pq{A.) (which is defined as the limit of (2.12)) 

n lyr^TATWj 2W-. 4t/^ + (A-A')^^°^ ^' (2_^3),3,, 



1 /"^ 

— / (ia; exp(ic<;A) 
27r J_oo 



Jo (a;) 



27r y_oo '2 cosh(;7a;) ' 

The ground state energy per site can be evaluated by using (2.9) 



^GS 

L 



U\-4 f dA Re - {A + i\U\)A Po{A) 



\TT\ I duj ex^{-\ulJ\) 
-1^1- / ,^ cosh(a;t/) 



(2.14) GSE 



Next we will determine excitations over the ground state (2.13) by means of the so-called 
"shift-function" method, which is explained in detail in [16] and [15]. 
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Charge-Triplet Excitation 

The lowest weight state of the charge triplet (77^ = —1) is is obtained by choosing 
M[ = ^ — 1. The allowed range (2.8) of the half-odd integers J'J^ is [— f , f ], so that there 



are f + 1 vacancies and thus two "holes" . We denote the spectral parameters corresponding 
to these holes by and Ag. The Lieb-Wu equations read 

r ~ ^ ^ ^"^^ A — A ^ A — A'* 

2L Re {arcsin(A, + = 27rJ^ + ^ ^(^^) - E ^(^M^) " ^2.15) BAECT 

Here the tilde indicates that the spectral parameters and the integers have changed as com- 
pared to their ground state distributions. Our convention is — = | for a. = 1 ... "I ^■ 
As compared to the ground state there is one more vacancy in the excited state and thus one 
more J^, which in our convention is the minimal allowed half-odd integer — ^. We denote 
the corresponding spectral parameter by A^-^ and note that in the thermodynamic limit 
^min ~^ —00. Subtracting (2.15) from (2.11) and using that A^ — A^ = ^(l) find 



(2.16) diffl 



We now define the "shift" function F^^j, for the charge-triplet according to 

PcAK) = "J^r , (2.17) F 

and take the thermodynamic limit of (2.16) 



p(A) y ]^/TZriATi\W J -00 4t/2 + (A-A' 

2 . ' (2-18) FO 



To order 0{L ^) p(A) is the same as Po(-'^)) that we can use (2.13) in (2.18), obtaining 
an integral equation for Fqj, 



^The results for energy, momentum, and scattering phase are of course independent of our choice of 
convention. 
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FcA^) = 1 - 



27r 



Equation (2.19) has the solution 

where 7(/i) is defined as'^^^ 



A- 



+ 7 



A-A^ 

2\U\ 



2\U\ 



(2.19) F2 



(2.20) F4 



7(/^) 



dco exp{iu>ii) 
Loo 1 + exp(|a;|) 

The excitation energy follows from (2.9) 



i In 



r(i±^ 



)r(i-f: 



r(i^)r(i + f) 



(2.21) fadtakh 



E. 



CT 



f 

J — c 



r-oo 

4 / dA Re 

-oo 

2 

+ ' 



A + ^|^7| 



Vl-(A + i|[/|)2 



F^^(A) 



where 



4^ [Re{0^^^A[7^} - 
ec«,(A?)+e,,(A^) , 

(\\ o du Ji{uj)cos{ujX) 



(2.22) ECT 



/q a; cosh(cj^) 
The momentum of the excitation is computed by means of (2.10) 

dA 2Re <^ " ^ i^CT(A) - X] {arcsin(Aj + i\U\)} - tt 



(2.23) es 



CT 



J — c 



Vl-(A + i|C/|)^ 



PcJA^)+p,t(A^ 



2; ) 



where 



PC 

pLW = - / 

Jo 



duj Jo(c<;)sin(a;A) 



(2.24) PCT 

(2.25) ps 



a; cosh (cut/) 

The quantum numbers of this type of excitation are (due to the lowest weight theorem) 
S — T]^ = —1, rj = 1, and the state is lowest weight of SO (A). The complete multiplet 
is obtained by acting with rj^ and {rj^Y, which leaves the energy but not the momentum 
invariant (77^ does not commute with the momentum operator). 

Charge-Singlet Excitation 



The charge singlet is obtained by choosing M[ — ^ — 2 and = 1. The allowed 

[.-2 L-: 
4 ' 4 



range of vacancies for the integers J'^ is [—^^, ^^]: so that there are again two holes with 
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corresponding spectral parameters Aj^,A2 in the distribution. The integer J'^ can take the 
value only. We denote the spectral parameter corresponding to the A;-A-string of length 2 
by K. The Lieb-Wu equations (1.4) for this type of excitation are 



2L Re {arcsin(A„ + = 2^3^ + ^^(%^) ' E ^(^^|^) 



2\u\ ' ^ ' m\ 

+ ^^^( "^1^ ^ ) (2-26) BAECS 



2L Re{arcsin(«: + 2^|t/|)} = ^^,,(^1^) - ^^^^(^^ 



In the thermodynamic limit they turn into integral equations for an F-function and a 
density p defined like in (2.17) and (2.12) 



/CO Lc \ ^ 

-oo 1^1 ^ 



=1 2|[/| ^ 

^ (2.27) BAECS2 

A. lA„ ,«:-A^ 



The second equation in (2.27) determines k as a function of A^. The zeroth order part of 
this equation is fulfilled identically, and the 0{L~^) part leads to the requirement 

K = — . (2.28) kappa 

The equation for F^g can then be solved by Fourier transforming 



where 7(/u) is defined in (2.21). Energy and momentum can be determined by using (2.9) 
and (2.10). We find 

Ecs = + e.M) , ^p^^ 

Pes = PLi^i) + PLi^2) , 

where e^^{\) and Pcwi^) given by (2.23) and (2.25), and where = t^+PcwW- Thus 

we see that this excitation has the same energy as the charge triplet excitation considered 
above, whereas the momentum differs by tt. Charge and spin quantum numbers are rj = ~ 
S. 

Spin-Charge Scattering States 
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We now choose M[ — f — 1 and Ai^ — N^ — 2M' — 1. is the number of "elementary" 
kj^s in a solution of (2.4)-(2.6). The number of vacancies (2.8) for the integers J'^ is ^, so 
that there is one hole with corresponding spectral parameter A'' in the distribution. The 
spectral parameter of the one real k is denoted by k. The Lieb-Wu equations read 



L 



2L Re {arcsin(A„ + i\U\)} = 2nJ-+± ^(^^W^) " ^(^TT ) + ^(^1^)' 

13=1 II II II 

(2.31) BAEes 

The second equation in (2.31) determines A; as a function of the integer /, which has range 
[— -^Y^, ^^yi]. In the thermodynamic limit the first part of (2.31) turns into an integral 
equation for an F-function 

2.F,M ^ - £ .A' W + "l^^) - "(^j , (2.32) 

which has the solution 

^ /.X 1 f , A-sin(A:)A 1 1 /A-A^\ 

F,.(A) = - arctan j - j - jjl^ j . (2.33) Fes 

The excitation energy of the charge-spin scattering state is given by 

A + «|f/| 



E^, = -2 cos(A;) + 4 y dK F^,(A)Re | -j= 



(A + z|[/|)^ 



+ 4Re{Vl-(A'^ + «|t^|)'} -2| 
= e,^(A'^) + e,^(A;) , 



(2.34) EES 



where ecjul-^) given by (2.23) and 

/,x ^ /,x ^ /"°° (io; Ji (a;)cos(a;sin(/c)) , , 

A; =2t/ -2cos A; +2 / ^^-^-f— ^ exp -^[/ . 2.35 ec 

Jo <^ cosh(a;L') 



The momentum is found to be 



n.=PL(A'')+P.«,(fc) , (2.36) PES 



where j^i}^ is given by (2.25) and 

/"°° <^<^ •^(<^)sin(a;sin(/^;)) , , 
"'-W^*-/ ^ coshU (2.37) PC 
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The quantum numbers of this excitation are —S^ — \ — S and —vj^ — \ — r] and the state 
is thus the lowest weight state of an S0{4) multiplet of dimension 4. The other states of the 
multiplet are obtained by acting with , r]\ and both. 

Spin-Singlet Excitation 

The spin-singlet excitation is characterized by taking M[ — ^ — 1, M.^ — 2, and 
= 1. The total number of electrons is A''^ = L, so that there is no charge associated with 
this excitation. The Lieb-Wu equations are 

2L Re {arcsin(A, + ^\U\)} ^ 2nJ-+ ^(%^) + E ^(^7^)' 

13=1 ' ' i=l ' I 

L k, = 2.1 + T O r^';^-^" ) + 9(5^^%^) (2.38) BAESS 

13=1 ' ' ' ' 



AC — sin(fc) 



^ \U\ 

j=i I I 



) = 27rJ^ . 



Here k is the spectral parameter corresponding to = 1, and k^,k2 are the momenta 
of the real /c's. There arc f — 1 vacancies for the half-odd integers J'^, so that there arc 
no holes in the distribution. The integer is found to be by (2.8), which implies that 
K = ^Hki)+sin{k2) ^Yie third equation of (2.38). The first equation in (2.38) once again 
turns into an integral equation for an F-function in the thermodynamic limit 



/ N 1 \^ / , A-sin(A;.)A 
^ssi^) = - 2^arctan (^exp(7r ) J 



1 A / A - sin(A;,) \ (2.39) FSS 

2^ - 



Energy and momentum are given by 

Ess = ^swik-i) + ^swik-i) , 
PSS ^ Pswikl) + Pswih) , 



(2.40) EPSS 



where es^(fc) and Pswi^) given by (2.35) and (2.37). The quantum numbers are S — — rj, 
so that the excitation is a singlet of 5*0(4). 

Spin-Triplet Excitation 

The last type of excitation in the two-particle sector is the spin-triplet excitation. The 
quantum numbers of the lowest weight state are r] — 0, — —1, S — 1. It is obtained 
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by choosing M[ — ^ — 1 and M.^ — 2, which leads to a total number of electrons — ^. 
There are ■§ — 1 vacancies for the half-odd integers J'^, so that there are no holes in the 
distribution. The Lieb-Wu equations are 

2L Re {arcsin(A„ + ^\U\)] = 2nJ'^+ ^(^^) + E ^( |^| )^ 

(2.41) BAEST 

/3=1 ' ' 

The F-function is found to be identical to ^^^(A), which implies that energy and momentum 
of the spin-triplet and spin-singlet are the same. 



3. QUASIPARTICLE INTERPRETATION FOR U < 

All of the excited states construced above have a natural interpretation as scattering 
states of four elementary excitations ("quasiparticles"), which form the fundamental repre- 
sentation of SU{2)gp^^ X S'C/(2)gj„, i.e. they are grouped into doublets of the spin and r]-S\J{2) 
respectively. The chargeless spin-carriers ( "spin- waves" ) have the dispersion 



Jo 



doj Jo(cc;)sin(cj sin(A;))e-'^l^l 



uj cosh(a;f/) ,„ ^ . 

irri 3.1) aepc 

e,^{k)^2U -2cos{k) + 2 ^ . 

Their quantum numbers are rj ~ 0, S — 1,5^ = ±|, which corresponds to the representation 

(i,0) of SU{2),^,^xSU{2),,,. 

The spinless charge-carriers ( "charge- waves" ) (one particle, one hole) have the dispersions^^^l 

(W ^ r ^o(^)sm(^A) ^.hr.. 
Jo 



ec«,(A)=2 



oo 



CO cosh(a;?7) ,^ ^, 

^ T r \ r w (3-2) aeps 

dcu Ji[u!)cos{ujX) 



00 cosh.(uU) 



Their quantum numbers are 8 = 0,7]=-^, r}^ = ±| and they form the (0, |) representation 
of SU{2),^^ X SU{2),,^. 

Like in the case of the spin-| Heisenberg antiferromagnet the quasiparticles are "confined", 
i.e. they do not exist as one-particle excitations. This is closely related to the fact that the 
symmetry of the Hubbard hamiltonian is 50(4) and all half-odd integer representations of 
SU{2) X SU{2) (and thus the fundamental representation formed by the four elementary 
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excitations) are not present^. The two-particle scattering states of the four quasiparticles 
arc easily seen to reproduce exactly the twelve excitations of section 2. 

(i) Scattering of two spin-waves 

The dispersion relations for the scattering states of two spin- waves follow from (3.1) to 
he E = e^,^(A;J + e^^ik^), P = p,^{k^) +Psv,ih)^ where k-^ and k^ are the rapidities of 
the spin-waves. The S'0(4)-representation is simply (|,0) ® (^,0) = (1,0) © (0,0). By 
inspection we see that (1,0) is identical to the spin-triplet excitation of section 2, and 
(0, 0) to the spin-singlet. 

(ii) Scattering of two charge-waves 

Scattering of two charge-waves leads to the 5'0(4)-representation (0, |)(8)(0, |) = (0, 1)© 
(0, 0). The energy of all four states follows from (3.2) to be £^ = e^^{Ai) + e^^{A2), where 
and A2 are the rapidities of the charge- waves. The lowest weight state r^^ = — 1 of 
the triplet (0, 1) contains two holes and has momentum P = p^^[K^) + p^^i^Kl^). It is 
identical to the charge-triplet state constructed in section 2. The other two states of 
the triplet have momenta P^.^^ = + P^cwi^2) and P^.^^ = p^^(AiO +P?^(A^') 

respectively, which is also in agreement with the results of section 2. The (0, 0) state is 
seen to be identical to the charge-singlet constructed above. 

(iii) Scattering of spin-waves and charge-waves 

The relevant S'0(4)-representation for this process is (|,0)(8)(0, 1) = (|,|). The lowest 
weight state of the multiplet has quantum numbers r]=^ = S, S^ = —^ = ri^ and 
a dispersion E = e^^(A'*) + e^^^k), P = p^^{A^) + Psw{k), where k and A'^ denote the 
rapidities of the spin-wave and charge-wave respectively. This state is clearly identical 
to the spin-charge scattering state of section 2. The other states of the multiplet are 
identical to the ones found in section 2 as well. 

A similar analysis can be carried out also in the sectors with four, six, etc. particles. All ex- 
citations have a natural interpretation in terms of scattering states of the four quasiparticles. 
The general proof of this statement is given in the appendix. 

4. Scattering Matrix 

In quantum mechanical scattering theory the scattering matrix can be extracted from 
the asymptotics of the wave-function of the scattering state^^^' . The boundary conditions 
of the quantum mechanical problem are free. This is in contrast to the periodic boundary 
conditions imposed in the Bethe Ansatz solution. In [23] it was shown how to modify 
quantum mechanical scattering theory to accomodate for this fact, and a general method 
for extracting the exact S- matrix from the Bethe Ansatz solution was given. In [20,21,24,25] 

^This is completely analogous with the spin-i Heisenberg model, which is 0(3)-symmctric for even chains, 
whereas the elementary excitations form the fundamental representation of SU{2)^'^^'^^^ . 
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this method was apphed to evaluate the exact S-matrix for Heisenberg models. Here we will 
apply this method to the case of the nested Bethe Ansatz and the Hubbard model. The 
S-matrix for parity-eigenstates in one dimension is simply a phase factor. The two-particle 
scattering-phase is equal to the phase obtained by moving particle 2 through the whole 
interval (one-dimensional box) in the presence of particle 1 minus the one-particle phase- 
shift, which is equal to the phase picked up by moving particle 2 through the box when 1 
is absent. The two-particle phases can be computed from the Lieb-Wu equations (1.4) by 
using our results for the two-particle excitations above. However there are no one-particle 
excitations present and thus it is impossible to evaluate the one-particle phase-shifts by 
means of the Bethe Ansatz. However it is possible to evaluate the relative phases between 
two-particle excitations by means of the Bethe Ansatz. Thus it is possible to evaluate the 
exact S-matrix up to an overall constant factor, which is the difference between a suitably 
chosen reference phase, and the true one-particle phase shift. In other words the method of 
[23] allows the evaluation of the logarithmic derivative of the S-matrix. In general this leads 
to the subtle issue of determining the overall phase of the S-matrix^^l , which for models like 
the spin- 1 Heisenberg antiferromagnet cannot be resolved in the framework of the Bethe 
Ansatz. For the case of the Hubbard model we are in much more fortunate position: In 
the limit C/ — > oo the charge sector of the repulsive Hubbard model reduces to scattering 
of free fermions. This can be seen by directly evaluating the many-body wave-functions in 
this limit'^^1. Thus the S-matrix (in the specific sector) must reduce the the S-matrix for 
free fermions (which is 1) in the limit U ^ oo, which allows us to fix the constant. In the 
U ^ oo limit of the attractive Hubbard model the charge-triplet must reduce to a scattering 
state of hard-core bosons with S-matrix —1. 

The complete S-matrix for the Hubbard model is 16 x 16 dimensional and blockdiagonal. 
It breaks up into 4 blocks (due to conservation of spin and charge quantum numbers), 
describing scattering of spin- waves on spin- waves, spin- waves on charge- waves, charge- waves 
on spin- waves, and charge- waves on charge- waves respectively 

fS^M \ 

s^M 
V smJ 

Charge-Charge Scattering 

The four parity eigenstates are the charge-singlet and -triplet states. As the scattering 
phase is the same for all members of the triplet we will only consider the lowest weight state. 
The general form for the S-matrix of quasiparticles in the fundamental representation of 
SU{2) is 



(4.1) SM 
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where a,b^c^d label the colliding quasiparticles, and A^, A2 are the corresponding spectral 
parameters. The phase factors for singlet and triplet are 

S, = S, + S, , S, = S,-S, . (4.3) ts 

We will now use the method of [23] to determine and S^. We start by considering the 
charge-singlet excitation. The phase for moving particle 2 through the ring in presence of 
particle 1 is 

2 A'* — A A'* — A'* 

S'M^i,^2) = Re {arcsin(A^ + m} + E^(^7n^) " ^(^TTt) 

f3=i ^1^1 ^1^1 (4.4) dl2cs 

\h _ \ h 

In order to obtain the true phase shift we now ought to subtract the phase ^chargei^-i) 
obtained by moving particle 2 through the ring in the absence of particle 1 (i.e. the one- 
particle phase-shift). Due to confinement of the quasiparticles it is however impossible to 
compute 5 charge- Instcad wc wiU subtract a reference phase (^0(^2)) which is chosen such that 
the phase ^^^(Ai , A^) — 5o(-^2) is a function of A2 — Aj only. The difference S^harge ~ ^0 leads 
to an overall phase- factor for the S- matrix (4.1), which will be fixed below. The reference 
phase is 

2 A'* — A 

5o(A^) = -2L Re {arcsin(A2^ + t\U\)} + ^^( , (4.5) dO 

which leads to the following phase-shift for the charge-singlet 



^C5(Ai,A2)=^cs(A^A2)-5o(A' 



2 

41 

dA 



/oo 

where /i^^ = 2\u\ ^ ^ 0' ^i-^) ~ 2arctan(a;), and 7(/i) is defined in (2.21). 

In order to determine 5*^ we have to consider the charge triplet excitation. The phase- 
shift for moving particle 2 through the ring in presence of particle 1 is 

L 

2 A'* — A A'* — A'* 

SM^l^2) = -2^ Re {arcsin(A^ + ^|t/|)} + ^(^Tn^) ' ^(^TTt) 

p=i ^1^1 ^1^1 (4.7) dl2ct 

A'* - A 
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The extra tt is due to the fact that the J'^ are half-odd integers for the charge-triplet but 
integers for the charge-singlet and the ground state, which will be used for reference purposes. 
The phase-shift is found to be 



(A^A^) = 55V(A^A2')-5o(A^) 

hk-h^ (4.8) dct 

= 27r Fct{^2) mod 27r = tt - mod 27r , ^ ^ > . 

Using (4.8) and (4.6) we can now determine the phase-factors 5'^ and 



Ss = exp(z5c5) 



f^4 + i r(i^)r(i + ^) 

r(i±i^)r(-\ 



2 

_ r(^^)r(i - ^) 



(4.9) sec 



r(i^)r(i + ^) ■ 



Thus the final result for the 4x4 block of the complete S-matrix describing scattering of 
Charge-waves on Charge-waves is found to be 



SAk\,K\) = - '\ ^ ^ ^ ^ (^^7+— nl . (4.10) Sec 

""^ ' ^' r(i^)r(i + ^) V/«4 + ^ + i ) 

where I and H are the 4x4 identity and permutation matrices respectively. This S-matrix 
is as a function of spectral parameter (but not of momentum) identical to the well-known S- 
matrix for the spin-| Heisenberg antiferromagnet'^^'^^l and the SU{2)i Wess-Zumino-Witten- 
Novikov (WZWN) model^^^J (this S-matrix also describes scattering of particles in the SU{2) 
Gross- Neveu model^^^', and spinous in the Kondo model ^^^'^^1 ) . 

Scattering of Charge on Spin and Spin on Charge 

The relevant excitation for these processes is the spin-charge scattering state of section 
2. We first consider the charge-wave to be the active scatterer. The phase for moving the 
charge-wave with spectral parameter A'^ around the ring in presence of the spin-wave with 
spectral parameter k is 

5l\{K\k) = -2L Re{arcsin(A^ + ^|^|)} + ^g(^^lz^) + ~ 'i""^^^ )- (4-11) dl2es 

/3=1 I I II 

The reference phase is once again given by (4.5), which yields the following phase-shift for 
charge-spin scattering 

5..(A^ k) = 5f,{K\ k) - 5,{K^) = 27r F^,(A^) mod 27r 

, , TT A''-sin(A;) ^ (4.12) des 

= 2 arctan(exp(7r;U3)) - - , /X3 = > . 
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The corresponding 4x4 block of the S-matrix is 

(^. ,^^,1+1^2^, (4.13) Ses 

Now we consider the spin-wave as the active scatterer. The total two-particle phase shift is 



To get the phase-shift we ought to subtract the one-particle phase-shift for a spin-wave 
Sgpin{k), which is unknown due to confinement. By the same method as above we can, 
however, determine the phase-shift up to a constant. The reference phase is Sglk) — Lk — 
Z]/I=i ^(~~^^]~^)) which yields the following result for 

5,^iA^k)^Sf^iA\k)-~S,ik) 

~ J-J C/2 + (sm(A;)-A')2 ''^^^ ^ \U\ ^ (4.15) dse 

, , TT sin(/c) — A'^ 
= 2 arctan(exp(7r//2)) - , 1^2^ ^y^^^ > • 

This is a priori correct up to a constant, which we will now show to be zero. We know that 
5^g{sm{x),x) = 5g^{sm{x),x), as in this case the quasiparticles are at rest relative to each 
other. Inspection of (4.15) and (4.12) now shows that the constant must vanish. As a result 
we see that it does not matter whether spin or charge is taken to be active. The form of the 
S-matrix (4.13) is the same, only the definition of /i changes. 

Spin-Spin Scattering 

In the spin-spin sector the S-matrix is once again of the form (4.2). We first determine 
the phase-shift for the triplet, as the J^-'s are integers both for the spin triplet and the spin- 
charge scattering state we used to determine the reference phase in the spin-sector. The 
phase for moving particle 2 through the ring in presence of particle 1 is 

Sf^k,, k,) = Lk,-J2 . (4.16) dl2st 



a=l 



Subtracting the one-particle phase-shift for spin- waves 5q(A;2) we arrive at the following 
expression for the triplet scattering phase 



sin(A;2) — sin(A;j) 



^sAh, = iM mod 27r , = — ' ,,,, > . (4.17) dst 
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The spin-singlet phase-shift is found to be 

^ssi^i-i ^2) — 7(a*i) ~ 2 arctan(yUi) -|- tt mod 27r , (4-18) dss 

which imphes the following form for the S-matrix in the spin-spin sector 

S„(sin(.f ), sin(.;)) ^ - ^n) . (4.19) Sss 

Last but not least we have to fix the overall factor of the S-matrix by considering the 
U — > —00 limit. In this limit the wave-function for the charge-triplet state reduces to the 
one for hard-core bosons. Therefore the charge-triplet phase-shift has to become tt in this 
limit (the S-matrix for hard-core bosons is —1, the one for free fermions 1). This fixes the 
constant factor. We see that our expression for the charge-triplet phase-shift (4.8) has the 
correct asymptotic behaviour, so that our constant factor is already the correct one. This 
completes our computation of the exact S-matrix for the attractive Hubbard model. 
The analytic structure of the S-matrix is related to the existence of bound states of the 
elementary excitations. The physical strips for spin-waves and charge-waves are defined by 
|lm(sin(A;))| < 2\U\ and |Im(A'^)| < \U\ respectively. Using these ranges in the expressions 
for the S-matrices S"^^, S^^ and S^^ we find poles at sin(A;2) — sin(/i;^) = 2i\U\ in the spin-spin 
sector, at A'' — sin(fc) = ±i|?7|^ in the charge-spin sector, and at A2 — A^^ = — 2i|[/| in the 
charge-charge sector. The pole in the charge-charge sector is at the boundary of the physical 
strip and corresponds to zero total energy and momentum. Thus it docs not correspond to a 
physical bound state. The analysis in the spin-spin and charge-spin sectors is more difficult. 
For large values of |C/| one can show that the poles either he on the unphysical sheet (and 
thus lead to anti-bound states, which due to non-normalizability drop out of the physical 
Hilbert space), or lead to excitations with fixed momentum (and thus cannot be interpreted 
as particles) and higher energies than the scattering states. We conclude that there are no 
bound states and the set of four elementary excitations is complete. 

5. Repulsive Hubbard Model 

The analysis of the Lieb-Wu equations (1.4) for the repulsive case is similar to the 
attractive case above. The analog of the logarithmic equations (2.4) - (2.6) was derived by 
M. Takahashi in [17]. The main difference is the form of the /c-A-strings: in the repulsive 
case A;^ = TT — arcsin(A'™ -|- i mil), k'^ = tt — arcsin(A'™ — i mil), whereas all other /c^'s 
are the same as in (2.2)(sce [17]). Equations (2.5) and (2.6) are the same for the repulsive 
case if we take \U\ ^ U. The analog of equation (2.4) for the case of repulsion is obtained 

^The sign depends on whether we consider the spin-wave or charge-wave to be the active scatterer. 
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by taking \U\ — > —U in (2.4). The boundaries of integers (2.8) are the same for [/ > and 
[/ < 0^. Energy and momentum are 



N^-2M' oo M: 



cos(/^J+4 J] J^Rel^l- (A'a" + «|C^|n)2} -2[/7V, + C/L , (5.1) rE 



1=1 n=l a=l 

N^-2M' oo 



P= ^ A;^. -^^(2Re{arcsin(A;^ + i|;7|n)} + 7r(n- 1)) . (5.2) rP 

j=l n=l a=l 

The half-filled repulsive ground state is characterized by taking A''^ = L, = ^, and 
then filling both the integers and the half-odd integers Ij symmetrically around zero. 
It is described by a set of two coupled integral equations and was found by Lieb and Wu 
in [5]. The two filled Fermi seas of the J^'s and the J^'s correspond to spin and charge 
degrees of freedom respectively. The excitation spectrum over this ground state has been 
previously studied by many authors'^^'^'^^'^'^^'^^'^^l^, although the Lie-algebraic structure of 
the excitation spectrum was not fully revealed. Like in the attractive case there are 12 
excitations in the two-particle sector. Unlike for the attractive case now two F-functions 
occur as the ground state consists of two Fermi seas. In order to explain our notation we 
will discuss the charge-triplet excitation in some detail and then quote the results for the 
other excitations. The Lieb-Wu equations for the repulsive ground state are 



a=l 

1=1 13=1 



(5.3) GSREP 



where Ii = I — -^y^, / = 1 . . . L, and = a — a = 1 . . . ^. In the thermodynamic limit 
(5.3) turn into the well-know set of coupled integral equations'^] 



(5.4) GSREP2 



11 /"°° 2U 

The ground state energy per site for the repulsive model is'^l 

= -[/-/ V/ J U^)JM ■ 5.5 rGSE 

L J_^ u) cosh(a;(7) 



^This is important for the validity of the completeness proof of [18] for both attraction and repulsion. 
^For a collection of reprints see [32] . 
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The charge-triplet excitation over the ground state is described by the following set of equa- 
tions {N^^L- 2, Ml = ^) 



L-2 

2 



a=l 



U 



(5.6) rBAECT 



1=1 



13=1 



2U 



There are two holes corresponding to spectral parameters sin(A;^), sin(A;2) in the distribution 
of the integers Ij. The half-odd integers are distributed symmetrically between and 
Our convention is Ij —Ij = — | and J\ — J\ = —\, where I- and J\ are the ground state 
distributions of "integers". Subtracting (5.3) from (5.6) and then taking the thermodynamic 
limit we can derive coupled integral equations for the F-functions F*(A ) = a'^"~^a 



3 ' kj 



1 



2U 



+ (sin(/c) - A)2 
AU 



^ 4C/2 + (sin(fc) - A') 



1 A-sm(k'}) 



(5.7) rFCT 



These equations can once again be solved by Fourier transforming 



1 A - siniB) 
^ct(A) = -1 + - 2^ arctan(exp(7r ^^7^)) 



sin(A;) — sin(A;^) 
2U 



(5.8) rFCT2 



where we have used (2.21). Energy and momentum of the excitation are given by 

EcAklk^) = AU + 2f2cos{k!;) + 2 f dk sm{k)F^j.{k) = e^) + ^hik^) , 

j=i 

PcAkl, 4) = -J2^J+ r dk FSAk) = pM) + pM) , 
.7=1 -^-^ 

dio Jo{uj)sm{io sin(A;))e-'^l^l 



(5.9) rEPCT 



where 



Phik) = 2 ~ ^ 



ef^{k) = 2U + 2 cos(A;) + 2 



cosh (a; t/) 
duj J^{uj)cos{io sin(A;))e-'^l^l 



(5.10) repc 



UJ 



cosh(a;?7) 
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Here e^(/c) and Ph{k) are energy and momentum of a quasiparticle carrying charge +e and 
no spin. This quasiparticle, which transforms in the (0, |)-representation of SU{2)^p-^ x 
SU{2)^f^, is called holon. The charge-triplet state is a scattering state of two holons with 
rapidities sin(A;f^) and sm{k2)- 

The charge- singlet is characterized by iVg = L, = ^ — 1, M'^ = 1. It is a two- 
parameter excitation with rapidities sm{ki) and sin(A;2) and leads to the following set of 
F-functions 

1 ^ A — sinf k'^) 
FSsi^) = -1 + - E arctan(exp(7r ^^)) 

"'"^ (5 11) rFCS 

Energy and momentum are found to be 

Pcsiktk^)^pM)+Pahik^) , 



(5.12) rPCT 



where eg^i^{k) = ei^{k) and Pah{k) = — vr + Pi^{k) are energy and momentum of the other 
member of the (0, |)-doublet. It is connected to the holon by action of r]^ and is called 
antiholon. The charge-singlet is thus the scattering state of one holon and one antiholon. 



The spin-charge scattering state is obtained by taking — L — 1, — ^ — 1 and filling 
the seas of integers Ij and J^. There is one hole in both distributions with corresponding 
spectral parameters A'^ and sin(A;'*). The F-functions for this state are of the form 



T^sr^^ ^ . / / A - sin(/c'^) , , 1 1 /A - A'^ 
F^^^(A) = - arctan(exp(7r ^ — -)) - - - —7 ' 



TT ' 2U A 27r' \ 2U 

1 1 / / sin(A;) - A'*,, 1 /sin(A;) - sin(A;'^) \ 

The dispersion relation is again of quasiparticle form 

E,sik\A') = e,{k'^) + eM'') , 
Pnsik\A')=PH{k')+Ps{A'') , 



(5.13) rFes 



(5.14) rEetas 



where €g{A) and Ps{A) are energy and momentum a quasiparticle in the (|, 0)-representation 
of SU{2)gp^^ X 5'[/(2)gj„, which is called spinon 



,(A) = 2 / 
Jo 



2 Jo u; cosh{uU) ' 
duj J-i^{u)cos{ujX) 



CO cosh(a;C/) 
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Spin-triplet and singlet states are scattering states of two spinons. The lowest weight state 
of the triplet is obtained from the Bcthc Ansatz by taking N^ — L, = ^ — 1. There are 
two holes with spectral parameters in the distribution of the A^,. The F- functions 

are found to be 

,.,1 1 ' /A-Af 



2 27r ^ ' \ 2U 



(5.16) rFST 
1 A sin(fc) 
FIAk) = - 2^ arctan(exp(7r — — '-)). 

The spin-singlet state is characterized by A^g = L, M-^^ — ^ — 2, — 1. Again there are 
two holes, and the F-functions are 

FUA) ^-±yj{ !^-JlL + 1 arctan(^^^ , 

27r ^ ' \^ 2t/ y TT ^ U ' ' (5.17) rFSS 

Energy for both triplet and singlet states are of quasiparticle form 

i?(A?,A^) = e,(At) + e,(A^) , 



m?,A2)=P.(At)+P.(A^), 



(5.18) rEPspin 



where eg{X) and ^^(A) are given by (5.15). The two-particle S-matrix for the repulsive case 
can be computed by the same method as for the attractive case. The phase-shifts are found 
to be 

CT CT\ 2) I r(i^)r(i + ^)y ' 2u 



5(yg — 27r F^g{k2) + TT — TT — 2 arctan(//4) + i ha 



r(i^)r(i - ^ 
r(i^)r(i+l^) 



SsT = 271 F|r(A^) = TT - Z In I ^ ], \rr > 

' '^ )T(i I ^) y 2?7 



= 27r -F|5(A2) = 2 arctan(//i) - i In 



r(i^)r(i - 
r(i^)r(i+^ 



(5.19) rphasel 

(5.20) rphase2 



71 sinfA;''^) A'^ 

Sr,s = 27r F^^,(A;'^) + tt = 2 arctan(exp(7r/X3)) ^ 2 ' = ^^^7 > 



^sr, = 27r F;(A'^) = 2 arctan(exp(7r/.2)) - 77 ,/^2 = > 



TT A'* — sin(A; 

2 '^2" 2C7 

Here we have fixed the over all constant by considering the charge-triplet in the C/ — > 00 
limit, in which the wave functions reduce to the ones for free fermions. Thus the S-matrix 



(5.21) rphaseS 
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must reduce to 1 in this limit. The complete S- matrix is then again of form (4.1), where now 



r(i^)r(i + ^) V/^i + ^ + ^ 

^.,(A^fe)^-^ ; + ^^"P|"^^| I, (5.23) rSse 

^''^ ^ 1 - z exp(7r//2) 

The analytic structure of the repulsive S-matrix can easily be inferred from the structure of 
the attractive one. There are again no physical poles of the S-matrix in the complex plane. 



6. Repulsion versus Attraction 

Attractive and repulsive quasiparticle spectra are related by an interchange of spin 
and charge degrees of freedom. To see this we first make the substitution k ^ n — k 
in the attractive case. This renaming of spectral parameters has of course no physical 
consequences. The spinon energy (5.15) for > is the same as the energy (3.2) of the 
charge waves in the attractive case. The same holds for the holon energy (5.10) and the 
spin- wave energy (3.1). On the level of the quasiparticle momenta these equivalences hold 
only up to constants. The spinon momentum is related to the momentum of the charge 
waves for ?7 < like ^^(A) = PcwW + f — P^wi^) ~ f • The spin- wave momentum Pgyj{k) is 
almost the same as the holon/antiholon momenta for f7 > Pcwi^) = p^(A;) -|- 1 = Pahi^) ~ f ■ 
These differences in the quasiparticle momenta are due to the fact that the transformation 
cji ^ (— l)''Cj-|. that interchanges spin and charge degrees of freedom does not commute 
with the momentum operator but changes the total momentum by vr. This can also be 
easily seen on the level of the two SU{2) algebras: the spin-S'[/(2), which commutes with 
total momentum, is transformed into the r]-SU{2), which changes momentum by tt. 
The S-matrices for attraction and repulsion after interchange of spin and charge labels are 
seen to be identical as functions of the uniformizing parameter 



7. Discussion 

In this paper we derived the quasiparticle interpretation for the one-dimensional Hub- 
bard model and evaluated the exact S-matrix describing the scattering of these quasiparticles. 
The classification of elementary excitations into spinous (which carry spin but no charge) 
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and holons/antiholons (which carry charge but no spin) clearly reflects spin-and charge sep- 
aration in the one-dimensional Hubbard model^. The S-matrix is a SU{2) x S'^(2)-invariant 
solution of the Yang-Baxter equation, which establishes the one-dimensional Hubbard model 
as a factorizable S'0(4)-scattering theory. Our S-matrix is very similar to the S-matrix for 
the SU{2) X SU{2) principal chiral model with Wess-Zumino term of level I'^^l The charge- 
charge and charge-spin sectors are also identical to the S-matrix for electron-electron and 
electron-kink scattering in the "exactly screened case" of the Kondo problem '^^1. Apart 
from its relevance in the field-theory limit the S-matrix is an important tool in the study of 
transport properties'^^'^^l , and plays a crucial role in the general theory of integrable models, 
where S-matrix and R-matrix are often proportional. In [37] a R-matrix for an embedding 
of the Hubbard hamiltonian in a commuting family of transfer matrices was derived. This 
R-matrix, which belongs to a fundamental model (where L-operator and R-matrix are essen- 
tially identical^^^^), is different from our S-matrix. It is quite natural to search for a different 
embedding of the Hubbard hamiltonian, using an R-matrix proportional to our S-matrix. 
The L-operator for such a construction would have to be different from the R-matrix. 
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The analogous statement for spin-waves and charge-waves in the attractive case is of course also true. 
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In this appendix we use the Thermodynamic Bethe Ansatz (TBA)!^^!'' to prove the 
vahdity of the quasiparticle interpretation (obtained by analyzing the two-particle sector) 
in the 2Ar-quasiparticle sector. We will discuss the U < case explicitly, the repulsive case 
can be treated in an analogous manner. The TBA-part of the analysis in the attractive case 
discussed here is very similar to the one for repulsion, which is presented in great detail in 
[17]. We refer to that paper for explanations and notations and merely give the results here. 
Starting point for the TBA are the Bethe equations (2.4) - (2.6). In the thermodynamic 
limit these equations turn into coupled integral equations for the particle-densities p, o"^, 
o"', and holc-dcnsitics p, a^, {n = 1, . . . cxo) describing the distribution of elementary fc's, 
A-strings of length n, and k — A-strings of length n respectively 



/, s _/, s 1 cos(A;) ^ 



dk 



2n\U\ 



n=l -°° 

- ^A^ - ^ r rlh 2n|^| 

""^^ > ~ 27r n2[/2 + (A - sin kf 



a;(A) = -Re 



1 



TT \^l-{K + in\U\f 



n2C/2 + (sin A; - A)^ 

oo 
m=l 

m 2n\U\p{k) 
^ n'^U^ + (sin(/c) - A)^ 



(^n(A) + <(A)) 



m=l 



(I.l) holes 



where 



^nm * f 



1 d 
2ndX 



dx'e^ 



A- A', 

• \U\ ' 



/(A') 



More precisely, (I.l) determine the holc-dcnsitics in terms of the particle densities. The 
particle densities (at finite temperature) are determined by minimizing the thermodynamic 
potential Q — E — TS — p-N [p is the chemical potential^, T the temperature, S the entropy, 
and E the energy following from (2.9)). This leads to the following equations for the ratios 



^For a review see [39] . 

^We will set p = 2U as we are interested in the half-filled case. 
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of densities C,{k) 

Hm) -- 

ln(l + r/JA)) = 



and 



-2 cos(A;) ~ A* _|_ 1 



n=l 

*in(i + — ; 



27r 



2n|[/| 



n2[/2 + (sin(A;) - A)^ 



ln( 



1 + 



Ml + ry^A)) 



n=l 

1 r 2n|[/|cos(fc) 
27r n2[/2 + (sin(A;) - A)2 

-4Re I yi- (A + m|[/|)2| - 2n// 



ln(l 



(1.2) TBA 



+ 



27r 



T 

dk 



n=l 



2n|f7| cos(A;) 



n2C/2 + (sin (A;) - A): 



The ratios of densities are related to the dressed energies K{k) (for elementary A;'s), e„(A) 
(for A-strings of length n) and e^(A) (for k — A-strings of length n) via C{k) = exp(^^), 



exp 



:^), and v'nW 



energies except e'i(A) are greater or equal to zero, whereas e\{K) < on the whole real line. 
This imphes that the attractive ground state is obtained by fiUing all vacancies for k — A- 
strings of length 1, which is how we constructed it in section 2. In the T = limit it is also 
possible to completely solve the system (1.2) for the dressed energies. We find 



■<{A)- 



In the zero-temperature limit all dressed 



e'i(A)=-2 r 

Jo 



duj Ji{u)cos{uA) 
cosh.{ujU) 



K{k) = 2\U\ -2 cos(A;) + 2 

lo ^ 

e;(A) = 0Vn>2 , e„(A) = V n > 1 



du J^(^')c()s(^' sm(A-))r""'l'^l 



(1.3) dresseden 



cosh{u>U) 



This is in perfect agreement with our analysis of section 2: making a hole at spectral pa- 
rameter A'* in the distribution of cr^(A) costs energy —e[{A^). This obviously corresponds 
to a charge- wave, and we identify €[ — — e^^. Introducing one elementary k increases the 
energy by K,{k), which is seen to be identical to the spin- wave energy egyj{k). The additional 
information we get out of (1.3) is that only these two dressed energies are nonvanishing. 
Thus the energy of any excited state is the sum over — e'^'s and «;'s. A similar analysis can 
be carried out for the total momentum of excitations at T = 0. Using the expression 




M' 



n=l a=l 



n=l a=l 



n=l 



for the total momentum it can be shown that only k — A-strings of length 1 and elementary 
/c's contribute dynamically {k — A-strings of lengths greater than 1 contribute 7r(n — 1)). 
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Thus we conclude that both energy and momentum of any excitation breaks up into sums 
over the quasiparticle energies and momenta. What remains to be shown is now that the 
quasiparticle interpretation (QPI) predicts the correct number of 5'0(4)-representations with 
spin S and ?7-spin rj in the sector with 2N quasiparticles (recall that there always are an even 
number). The prediction of the QPI for the number of lowest-weight scattering states of M.^ 
spin-waves with spin S (note that there is a constraint that ^ — 5* is always an integer) is 
simply 



Xi 



Ale 
2 



Ale 
2 



S-l 



This number has to be multiphed by the number of lowest- weight scattering states of 2N—M.^ 
charge- waves, which is 



X2 



2N-M, 



2N-M^ 



N- 



— r] — 1 



Are these numbers reproduced by the Bethe Ansatz ? The answer is of course "Yes". 
Excitations over the ground state are construced based on the allowed ranges of integers 
(2.8). By analyzing these equations we can determine the exact number of Bethe Ansatz 
excitations (i.e. lowest weight states of >S'0(4) multiplets). We fix the number of holes in the 
A'^ distribution to be 2A^ — M.^ and the number of elementary A;'s to be Ad^. The dispersion 
of such states is by the above arguments identical to the dispersion of a scattering state of 
Me spin-waves and 2N — M.^ charge-waves. If we further fix the values S of spin and 77 of 
ry-spin, the total number of Bethe states with these characteristica is given by 



E 



n=l 



X 



E 



E 



n 

n=2 



m=l ^mn^m 



X 



2N -M^ 



mn 

M' 

n 



2) ML 



These sums can simply be looked up in [18] or [40] and we find that the two factors exactly 
coincide with Xi and X2- This establishes the validity of the QPI in the 2A^-quasiparticle 
sector. 
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